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Abstract. The possibility of realization of a superfluid state of bound electron-hole pairs (magnetoexcitons) with spatially 
separated components in a graphene double layer structure (two graphene layers separated by a dielectric layer) subjected by 
a strong perpendicular to the layers magnetic field is analyzed. We show that the superfluid state of magnetoexcitons may 
emerge only under certain imbalance of filling factors of the layers. The imbalance can be created by an electrostatic field 
(external gate voltage). The spectrum of elementary excitations is found and the dependence of the Berezinskii-Kosterlitz- 
Thouless transition temperature on the interlayer distance is obtained. The advantages of use graphene double layer systems 
instead of double quantum well GaAs heterostructures are discussed. 
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INTRODUCTION 

It is believed that excitons may demonstrate superfluid behavior. In bilayer systems superfluid excitons consisting of 
an electron from one layer and a hole from the other layer behave as superconductive ones. Indeed, having separate 
contacts in each layer one can use excitons for a nondissipative transmission of an electrical current from the source to 
the load (Fig. [TJ. The effect can be realized in double quantum wells in GaAs heterostuctures [1]. If such a system is 
subjected by a strong perpendicular to the layers magnetic field and the total filling factor of the Landau levels V = 1, 
electrons that occupy quantum states in the zeroth Landau level in one layer couple with holes (empty states in the 
zeroth Landau level) in the other layer. Such pairs, called magnetoexcitons, are stable one, in difference with optically 
excited indirect excitons in double quantum wells Q|3|]. Stable excitons can emerge also in bilayers made of n-type 
and p-type two-dimensional conductors flllt]. In the latter systems the nesting of the Fermi surfaces of electrons and 
holes is required for the Bardin-Cooper-Schrieffer (BCS) pairing of n and p carriers. 

The discovery of graphene Ja, |7[] has risen the idea of use double layer graphene systems for the realization of the 
exciton superconductivity. This question was already investigated in a number of papers 

Graphene can be considered as a semiconductor with zero band gap. Electron energy spectrum of graphene contains 
two Dirac points that separate the electron and the hole subband. In a bilayer structure the Fermi levels of the layers 
can be adjusted independently by the gate voltage. If the Fermi level is in the conduction band of one layer and 
is in the valence band of the other layer we have a n-p bilayer. The electron-hole symmetry near the Dirac points 
ensures perfect nesting between the electron and the hole Fermi surfaces. Such a situation was considered in 181,191, 1 1011 . 
While the estimates in which the screening effects are neglected yield very high BCS temperature (hundreds of 
Kelvins), the screening may reduce the critical temperature down to 6 orders iflOfl . 

In ITll [l2Tl the superfluid transition in the rarefied gas of magnetoexcitons was considered. Rarefied means that 
the number of magnetoexcitons is much smaller than the number of states in the Landau level. Since the critical 
temperature in 2D is proportional to the concentration of the carriers, the rarefied gas of magnetoexcitons cannot 
demonstrate record critical parameters. The aim of this paper is to analyze the case of high magnetoexciton density 
in bilayer graphene systems. Here we do not consider the effect of screening. In quantum Hall systems screening is 
expected to be small due to finite gaps between Landau levels. 



3 




+\ 2 



A B 




FIGURE 1. Bilayer graphene system (a dielectric matrix (1) with two embedded graphene layers (2) inside a capacitor (3)) that 
transmits the current by magnetoexcitons. 

INTERLAYER PHASE COHERENCE IN A DOUBLE LAYER GRAPHENE SYSTEM 



In conventional quantum Hall bilayers the case of high magnetoexciton density is realized at zero imbalance of filling 
factors of the layers (Vi = V% = \ /2). The system at zero or moderate imbalance of filling factors can be regarded as 
an easy-plane quantum ferromagnet []13|. Here we extend the approach of Ref. [13] to a graphene bilayer system. 
The Landau level energies and eigenfunctions for graphene are given by the equation 11411 
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Here Vf s» 10 6 m/s is the Fermi velocity (for the Dirac spectrum vp is the matter parameter independent of the Fermi 
energy), i = y/fic/eB is the magnetic length, a = ± 1 is the valley index (that corresponds to the carriers with the 
pseudo-momentum close to the K = {2n /3a, 2k /3^/3a) and K' = (2n/3a,—2n/3V3a) Dirac points, correspond- 
ingly), the spinor index A (B) corresponds to the A(B) sublattice of graphene. The solution of Eq. (Q~|i yields the energies 
of the Landau levels in graphene 

E = 0, E ±N = ± T ^fV2N, (2) 



and the corresponding eigenfunctions ^N.a 
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Here X — k£ 2 is the guiding center of the orbit, Hn(x) is the Hermite polynomial, and L y is the size of the system in y 
direction. Approximation (1) is valid for \En\ <C 3f, where 3t w 8 eV is the half-width of the energy band of graphene 
(in the tight-binding approximation with the nearest-neighbor hopping energy t). We do not take into account the 
Zeeman splitting because it is much smaller than the energy distance between the Landau levels. 

Due to the spin and valley degeneracy the total number of quantum states in the Landau level is 4S/2n£ 2 (S is the 
area of the system). We define the filling factor of the Landau level as v = 2n£ 2 rif, where n/ is the concentration of 
occupied states. A completely filled level in graphene has the filling factor v = 4. In undoped graphene the negative 
levels are completely filled, the positive levels are empty and the zeroth level has the filling factor v — 2. 



The electron-hole pairing occurs between the carriers that belong to a partially filled Landau level (we will call it 
the active level). In particular, the pairing of electrons and holes belonging to the zeroth Landau level may take place. 
The many-body wave function that describes the state with such a pairing can be written analogously to 11311 
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Here /3 = (a, a) are the sets of valley and spin quantum numbers (below we will notate them by digits /3 = 1,2,3,4), 
a tkB ' s ^ e °P era t° r of creation of an electron in the z'-th layer in the active Landau level, <pp is the phase of the order 
parameter for the electron-hole pairing A^ = ( l P|a ] H ^a2,j3 W = (1/2) sin Ope'^P , the parameter Op is connected with 
the filling factors of the component j3 by the relation Vim p — (1 ±cos 0p)/2. The vacuum state |0) is the state with 
empty zeroth (active) and higher Landau levels in both layers. 
Note that the function (0 can be presented in another form 
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where ht x ^ is the hole creation operator, and the vacuum state \vac) is the state, where the zeroth Landau level in the 
layer 2 is empty, while in the layer 1 this level is completely filled. From (0 it becomes clear that the function (0 is 
just the analog of the BCS wave function. The state ((5]) is usually called the state with a spontaneous interlayer phase 
coherence. 

We imply that the Coulomb energy E c = e 2 /e£ (s is the dielectric constant of the matrix in which the graphene 
layers are embedded) is much smaller than the energy distance between the active and the nearest passive (completely 
filled or completely empty) Landau level and take into account only the Coulomb interaction between electrons in the 
active level. Below we consider the case of N = active level. The Coulomb interaction has the form 

He = ^EEM^MqJfrHO. ( 7 ) 

where Vi/(q) = (2ne 2 /eq)exp(— qd\i — i'\) is the Fourier component of the Coulomb potential, d is the interlayer 
distance, and 

(8) 
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is the Fourier component of the electron density. 

We consider the bilayer system situated inside the capacitor that creates an electrostatic field normal to the graphene 
layers. Varying the electrostatic field one can change the imbalance of filling factors of the layers. The energy of the 
system in the state (0 reads as 
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where 

are the energy constants that describe the direct interlayer interaction, the exchange intralayer interaction and the 
exchange interlayer interaction (erfc(x) is the complementary error function). In (0 V is the external gate voltage 
caused by the capacitor, and vp = (Vj a — V2.p)/2 is the filling factor imbalance for the component j3 (|Vfi| < 1/2). In 
(0 the interaction of electrons with the positively charged background is included. The result (0 can also be obtained 
in the standard mean-field approach 111511 . 

The interaction constants satisfy the inequalities > ^\ and W — ^fo + fl\ > 0. One finds that at V = the 
minimum of (0 corresponds to Vi = V2 = 1 /2, V3 = V4 = — 1/2. At such vp all the order parameters A^ are equal to 
zero and there is no electron-hole pairing. If the gate voltage is in the range 



Jo-Ji<eV<2W-J + Ji, 
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the minimum is reached for 



~ - 1 - eV - W - 1 
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and A3 7^ 0. 
For 

2W + ^o~,A < eV <W-,/o+A (12) 

the energy minimum corresponds to 

1 - eV-3W 

Vl = V2 = V3 = 2' V * = 2(w- A)+Jl y (13) 

and A4 7^ 0. Thus we conclude that an imbalance of filling factors is required for the electron-hole pairing and only 
electrons and holes with one j3 are involved in the pairing. In the special cases eV — W and eV = 3W the imbalance of 
the active component has zero value, and the order parameter for the electron-hole pairing reaches the highest value. 

The difference between the case of a quantum Hall bilayer in GaAs with the total filling factor v = 1 and the case 
of the graphene bilayer system is the following. There is only one component in the first case and an imbalance of that 
component increases the direct interaction energy. In the second case positive imbalance of some components can be 
compensated be negative imbalance of the other components. The latter situation is similar to one that takes place in 
V = 2 quantum Hall bilayers II 1611 . 



CRITICAL TEMPERATURE 

The state with moving superfluid pairs is described by the order parameter with a spatially dependent phase. The 
many-body wave function of the state with a constant superfluid magnetoexciton current in x direction reads as 

W=n( cos y<x /j +^ sin y<x /J )|o> (14) 

(the only part that corresponds to the active component is displayed). The part of energy that depends on 0o and Q is 
equal to 
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(Jo(q) is the Bessel function). At small Q Eq. (TT~5T > is reduced to E m f = S(const + p s oQ 2 /2), where 



psO = sin 2 1 



I6nel 



-exp^jertc^ ^ + ^"7 



(17) 



is the mean-field value of the superfluid stiffness (we take into account that Q = V<p). At finite temperatures the 
excitations reduce this quantity. Therefore, the superfluid stiffness p s depends on temperature. 

Since the gas of electron-hole pairs in bilayers is a two-dimensional one, the transition into the superfluid state is 
the Berezinskii-Kosterlitz-Thouless (BKT) transition. The critical temperature of the BKT transition is given by the 
equation 

T c = ^Ps(T c ). (18) 

To obtain the dependence p s (T) one should find the spectrum of elementary excitations. Extending the approach 
fnl [l8ll to the general case of an arbitrary angle between the wave vector q and V<p we arrive at the following 
expression for the energy of excitations 



E(q) = y / e q , e (eq, e + 27 q , e sin 2 0o)+cos0ov q , e , (19) 
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(i is the unit vector in x direction). Having the spectrum of elementary excitations one can compute the free energy 

F=E mf + TY^ln(l-exp(-^-)) (25) 
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and the density of the superfluid current 
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On the other hand, at Q — ► the superfluid current is connected with the superfluid stiffness by the relation 

U = ^P-Q. (27) 
n 

Using (l24l)-([26l), we obtain the following expression for the superfluid stiffness 
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where Nb{E) — [exp(E/T) — 1] is the Bose distribution function. For the spectrum E(q) = Eo(q) + hq\ (where v 
is the superfluid velocity, and Eg(q) is the spectrum at v = 0) the temperature part of Eq. ( 1281 yields the standard 
expression for the normal density 11911 . The present situation differs from [ 19] because the dependence of the spectrum 
on the superfluid velocity v = hQ/M (M = h 2 sin 2 9q/87ZP s q£ 2 is the magnetic mass of the electron-hole pair) is more 
complicate. 

The nature of such a difference is the following. In certain sense, superfluidity of electron-hole pairs in the bilayer 
can be considered as a kind of a counterflow superfluidity l20ll . Indeed, due to the Coulomb interaction electrons 
from the top layer are coupled with holes from the bottom layer as well as holes from the top layer are coupled with 
electrons from bottom layer. In other words, one can say about two species of bosons of different polarization. In such 
a two-specie system only a counterflow motion of species is possible. In the limit of low mangetoexciton density the 
second specie can be considered as an inert background and the system behaves as a one-component gas of interacting 
Bose particles. But at zero imbalance of the active component we have two equivalent species that move in opposite 
directions with the same velocities. In particular, it results in that the last term in ( fT9l ) goes to zero at 9q = JC/2. The 
counterflow superfluidity is a special kind of the two-component superfluidity. In similarity with the spectrum ( fT9l , 
the spectrum of elementary excitations of two-component superfluid system demonstrates a complicate dependence 
on the superfluid velocities (see, for instance, 112 ill ). 

Let us describe the behavior of the spectrum ( fT9b at different d/£ and Q = 0. At d — > the spectrum reduces 
to the quadratic one E(q) = h 2 q 2 /2M (for small q). According to the Landau criterium of superfluidity it means the 
absence of superfluidity. At d sa I the roton-like minimum emerges at the dependence E(q). The depth of this minimum 
increases under the increase of d/l and at some critical d = d c the minimum touches the q axis. Such a behavior of 
the spectrum means that at d > d c the system becomes instable with respect to a formation of the charge density wave. 
We imply the superfluidity does not survive in the charge density wave state. At 9q — n/2 (that corresponds to zero 
imbalance of the active component) d c w \.\15t. One can expect that at d — > and d — > d c the superfluid stiffness is 
suppressed, that results in a lowering of the critical temperature. 




FIGURE 2. The dependence of the critical temperature on the interlayer distance (d = d/fy, 1 - coupling in the zeroth Landau 
level, 2- coupling in the N = ±1 Landau levels. The quantities 7rp s o/2 are shown by dashed lines. 

The critical temperature obtained from Eq. ( fT8l l for the case of zero imbalance (6q = n/2) is shown in Fig. [2] One 
can see that, indeed, at d — > and d — > d c the critical temperature goes to zero. We also find that at intermediate d the 
critical temperature can be evaluated with a good accuracy from the mean-field superfluid stiffness. Since the stiffness 
is proportional to sin 2 0q = 1 — 4v 2 , the highest critical temperatures can be reached at zero imbalance of the active 
component (or at the filling factors of the zeroth Landau level Vi = 3/2, V2 = 5/2 and Vi = 1 /2, V2 = 7/2). 

The approach presented can be easily generalized for the case where electrons in the +N Landau level are coupled 
with holes in the —N level. This situation can be realized at large external gate voltage eV > 2\f7hvp j I. To obtain 
the critical temperature in the latter case one should take into account the additional factor fy{qi) = [Zw(g 2 £ 2 /2) + 
Ln-i (<? 2 ^ 2 /2)]/2 in the electron density operator (O (Ln(x) are the Laguerre polynomials). Respectively, the factor 
fx(k) appears under the integrals in (TToT i and ( 124-b . and the factor f^(q£) should be added into the definition (1231 . Here 
we consider the case of the ±1 active level. The mean-field superfluid stiffness in this case reads as 

2, «* 

16ne£ 

At all d/l the superfluid stiffness for the ±1 Landau levels is smaller than one for the zeroth level (at the same 
imbalance of the active component). For the ±1 active levels the critical interlayer distance is d c ~ 0.549^. Using the 
procedure, described above, we compute the critical temperature. The result is shown in Fig. [2] 

DISCUSSION 

We conclude that for the gas of magnetoexcitons in a bilayer graphene system the maximum critical temperature of the 
superfluid transition can be achieved for the coupling of electrons and holes in the zeroth Landau level and at special 
values of the imbalance of filling factors of the zeroth level (Vi — V2 = ±1 and V\ — Vi = ±3). 

For the formation of magnetoexcitons the Coulomb energy should be much smaller than the energy distance between 
an active and the nearest passive level. The same condition allows to neglect the effect of screening on the critical 
temperature. Indeed, screening can be taken into account by a substitution a g-dependent dielectric function e(q) 
instead of the dielectric constant e into the Fourier components of the Coulomb potential (see, for instance, B22I0 . In 
the case considered here the difference e(q) — e is small by the parameter E c j (E\ — Eq). 

For the graphene system the inequality E c <C E\ — Eq is just the condition on the dielectric constant e 3> £ c ~ 1 -5. 
At large e this inequality is fulfilled at all values of magnetic fields (of course, it is implied that the temperature is 
much smaller than E\ — Eq). The situation differs significantly from one that takes place in quantum Hall bilayers in 
GaAs heterostructures. For the carriers with the quadratic dispersion the energy distance between the Landau levels is 
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FIGURE 3. The dependence of the critical temperature on the inverse magnetic length at fixed d. Solid curve - magnetoexcitons 
in the zeroth Landau level, dashed curve - magnetoexcitons in the iV = ±1 levels. 



proportional to B (co c = eB/m^c, where is the effective mass) and the condition E c < Ti(O c is fulfilled only at large 
magnetic fields (B > 10 T for GaAs). The graphene bilayers require much lower B. 

It is interesting to analyze how the critical temperature depends on B. Let we have a bilayer system with a given 
interlayer distance d. We vary magnetic field adjusting simultaneously the gate voltage (to keep the imbalance close 
to the optimal value) and try to achieve the highest critical temperature. Presenting the dependence of Fig.[2]in e 2 /ed 
units we obtain (Fig. [3) the dependence of T c on l^ 1 (or on B 1 / 1 ). One can see that the critical temperature is restricted 
from above by T c ,max ~ O.Ole 2 / ed, and the maximum is reached for t in the interval [d -j- 2d]. To achieve the critical 
temperature T c = 1 K we may take the bilayer with the interlayer distance d as large as 400 A (for e = 3.9 (SiCh)) and 
apply comparatively low magnetic fields B > 0.4 T. At the same time, our estimate shows that extremely high magnetic 
fields are required if we want to realize a superfluid state of magnetoexcitons at high temperatures. For instance, for 
d w 10 A we have T c max w 4QK, but the magnetic field should be B > 600 T. 

There is another reason that makes the case of low magnetic field important. Low fields mean large critical interlayer 
distances. As was shown in I23II24I1 . in the bilayer system the genuine superconductivity is reached in the limit of zero 
interlayer tunnelingQ If the systems with a finite tunneling is used for the transmission of the current from the source to 
the load (Fig. [TJ, the dissipation is nonzero. The dissipation is connected with that the state becomes nonstationary at 
nonzero difference of electrochemical potentials between the layers. This difference is required to provide the current 
in the load circuit. The power of losses is proportional to the square amplitude of the interlayer tunneling j23l 12411 . A 
partial solution of this problem was proposed in [26]. In the setup, considered in 12611 . a stationary state with moving 
electron-hole pairs can be realized, but, regrettable, the setup B26I1 (and any other setup based on the same idea [24]) 
cannot be used for the transmission the current from the source to the load. The amplitude of tunneling decreases 
exponentially with the increase of the interlayer distance, and for d w 400 A the effect of tunneling can be completely 
neglected. 
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